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ABSTRACT 

2 ■ The color-magnitude (CM) diagram of cool dwarfs and brown dwarfs based on 

O ■ the recent astrometry data is compared with the CM diagram transformed from 

^ ■ the theoretical evolutionary tracks via the unified cloudy models (UCMs) of L and 

D ■ T dwarfs. A reasonable agreement between the models and observations is shown 

^ ■ for the whole regime of ultracool dwarfs covering L and T dwarfs, and this is 

^ ■ achieved, for the first time, with the use of a single grid of self-consistent nongray 
model photospheres accommodating dust cloud (UCMs; 700 < T^s < 2600 K). 

^ ■ A distinct brightening at the J band in the early T dwarfs revealed by the recent 

^ . parallax measurements is explained as a natural consequence of the migration 

T— I . of the thin dust cloud to the inner region of the photosphere and should not 

^ . necessarily be evidence for Burgasser et al.'s proposition that the dust cloud 

. breaks up in the L/T dwarf transition. Also, the rapid bluing from the late L 

^ . to the early T dwarfs is a direct result of the transition of the thin dust cloud 

Q_|. from the optically thin (r < 1) to thick (r > 1) regimes while Lboi and Tcs lower 

Q . only slightly. Thus, the theoretical evolutionary models, the cloudy models of 

+i . the photospheres (UCMs), and the observed fundamental stellar parameters are 

ci . brought into a consistent picture of the newly defined L and T dwarfs. 

^ , Subject headings: molecular processes — stars: atmospheres — stars: fundamen- 

^ . tal parameters — stars: late-type — stars: low-mass, brown dwarfs — 



1. INTRODUCTION 

Since the discovery of the genuine brown dwarf G1229B (Nakajima et al. 1995), a large 
number of ultracool dwarfs have been discovered (e.g. Strauss et al. 1999; Kirkpatrick et 
al. 2000; Burgasser et al. 2002a). The newly discovered ultracool dwarfs are now classified 
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into the new spectral classes L and T (Geballe et al. 2002), which are characterized by 
the very red colors due to dust extinction and by the strong bands of methane and water, 
respectively. Extensive efforts have been made to interpret these new objects in terms of 
atmospheric chemistry including formation of dust clouds, as reviewed recently by Burrows 
et al. (2001). Although L and T dwarfs appear to be quite different, we proposed that they 
can be understood consistently with the unified cloudy models (UCMs) in which a thin dust 
cloud is formed always near the dust condensation temperature and hence will be located 
relatively deep in the photosphere for the cooler T dwarfs while will appear in the optically 
thin regime in the warmer L dwarfs (Tsuji 2002). The presence of the dust cloud deep in 
the photospheres was also shown by an apphcation of the planetary theory (Marley et al. 
2002). It is interesting that the presence of the thin dust cloud deep in the photospheres of 
ultracool dwarfs has been concluded from the quite different approaches. 

The nature of the cloud, however, is by no means clear yet. One major interest is if 
the cloud is subject to the meteorological activities familiar with the planets in the solar 
system. Such a possibility has so far been suggested by observing photometric variabilities 
which may be due to the inhomogencity of the dust clouds (e.g. Bailer- Jones & Mundt 2001; 
Martin, Zapatero Osorio, & Lehto 2001). Also, an interesting finding of the recent parallax 
measurements is that the J magnitude shows a large brightening in the early T dwarfs (Dahn 
et al. 2002). This result was interpreted as evidence for the disruption of the dust cloud in 
the L/T dwarf transition (Burgasser et al. 2002b). 

Before such a possibility is explored, however, it is important to remember a more 
classical application of the CM diagram; namely as a touchstone of stellar models. For 
this purpose, the observed CM diagram of ultracool dwarfs including brown dwarfs is now 
accurate enough to be confronted with the theoretical evolutionary models (e.g. Dahn et al. 
2002). On the other hand, theoretical evolutionary tracks of the substellar mass objects have 
been discussed by Hayashi & Nakano (1963) already in the 1960's, and recent developments 
with improved input physics have been discussed by Burrows et al. (1997) and by Chabrier 
et al. (2000). Given that the observed CM diagram and theoretical evolutionary models are 
both reasonably well established, a missing link between them is a realistic model photosphere 
to be used for the conversion of the fundamental stellar parameters to the observables. We 
show in this letter that the missing link might be found in the UCMs (Tsuji 2002). 

2. COLOR-MAGNITUDE DIAGRAM 

We applied the UCMs to convert the Mboi and T^s of the evolutionary models to more 
easily observable monochromatic absolute magnitudes and color indices; we discuss Mj and 
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J — X as an example. In the UCMs, we assume that the dust grains formed at the conden- 
sation temperature (Tcond) are in detailed balance with the ambient gaseous mixture so long 
as the grain sizes are smaller than the critical radius (r^) and hence can be sustained in the 
photosphere. Dust grains soon grow to be larger than at a slightly lower temperature 
which we refer to as the critical temperature (Tcr) and the dust grains are stabilized. Such 
stable grains may segregate from the gaseous mixture and can no longer be sustained in the 
photosphere. The exact value of is difficult to know but we assume it to be quite small 
in the sub-micron regime (e.g. ~ O.Ol/itm), since astronomical grains prevailing in the 
universe (and hence these grains are already larger than Tcr) are well above this size as can 
be inferred from the reddening law. For this reason, only dust grains smaller than survive 
in the temperature range of Ta ^ T < Tcond- This means a formation of a dust cloud with a 
rather high temperature (note that Tcond ~ 2000 K) independently of T^s. Since T ^ Tcfj at 
Tross ~ 1) the dust cloud appears in the optically thin region (tross < 1) in L dwarfs whose 
Teff's are relatively high and in the optically thick region (tross ^ 1) in T dwarfs whose 
Tefj's are lower. As a result, the thin dust cloud formed in the photosphere of ultracool 
dwarfs moves from the optically thin region in L dwarfs to the deeper optically thick region 
in T dwarfs. This migration of the dust cloud in the photospheres of ultracool dwarfs has 
observable effects to be reflected on the observed CM diagram, as we will show below. 

In discussing model photosphere, one important input parameter is the chemical com- 
position. After the first version of the UCMs (Tsuji 2002), an important revision of the solar 
carbon and oxygen abundances was proposed by AUende Prieto, Lambert, & Asplund (2001, 
2002), who showed that the C and O abundances are about 50% lower than those by Anders 
& Grevesse (1989). In view of the decisive role of C and O in the chemical equilibrium in 
cool dwarfs (Loddcrs & Fegley 2002), we have computed a new grid of the UCMs with the 
new C and O abundances. The other abundances are by Anders & Grevesse (1989) and we 
assumed that the abundant refractory elements Fe, Si, Mg, and Al arc condensed into iron, 
enstatite (MgSiOa), and corundum (AI2O3). The extinction coefficients of these dust species 
are evaluated with the use of the measured optical constants and they do not depend on the 
grain size so far as the grains are as small as we have assumed (Tsuji 2002). 

To convert Mboi predicted from the evolutionary models to Mj, we apply the bolometric 
correction for the J magnitude defined by 



BCj = Mboi - Mj, 



(1) 



with 
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Mj = -2. Slog 
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where is the emergent flux calculated from our UCM with the use of the molecular linelist, 
Ci and C2 are constants which depend on the zero-points of the magnitude systems, R is 
the stellar radius, and dio is the distance to the object corresponding to 10 pc. We used the 
filter response function given by Persson et al. (1998), which is normahzed as 



where C — Ci — C2. We first apply eqn.(5) to determine C by the use of the model flux 
of Vega with T^s = 9550 K, log g = 3.95, and Kiicro = 2 km s"^ (Kurucz 1993). The value 
of BCj for Vega is -0.22 from mtoi = -0.22 (Code et al. 1976) and J = 0.0. We found 
C = —7.72, with which eqn.(5) is apphed to the model fluxes of L and T dwarfs. 

The resulting values of BCj against T^s are shown in Fig. la for the four values of T^, 
the fully dusty models (case B; = ^surface), dust segregated models (case C; = TLond 
and hence effectively dust- free), and the blackbody (BB) radiation. We assumed logg' = 5.0 
which is the median value of log g extending from 4.5 to 5.5 in brown dwarfs (Fig. 9 of Burrows 
et al. 1997) . The changes of log ^ by ±0.5 at T^s = 1600 K {T^, = 1800 K) , for example, result 
in the changes of BCj by about ±0.15 (Fig. la). The BCj can be regarded as a generalized 
color index and the larger BC^ of the case C compared with the case of BB for the lower 
Teff is due to the large nongray opacity which is relatively transparent in the J band region 
compared to the longer wavelength region dominated by the increasingly larger absorption 
due to H2 CIA (collision-induced absorption) as well as CH4 and H2O ro-vibration bands. 
If the thin cloud is introduced, the dust extinction at the J band region first increases from 
the early L to late L dwarfs and BCj decreases. The dust column density in the observable 
photosphere attains the maximum value of 7.5 x 10~^ gram cm~^ at T^s ~ 1500 K which 
corresponds to the late L dwarfs. In T dwarfs, however, the cloud moves into the optically 
thick region of the photosphere, and the observable effect of the dust extinction decreases 
in the L/T transition objects and further in T dwarfs, resulting in the re-increase of SCj. 
The values of BCj depend on T^^, which is essentially a measure of the thickness of the dust 
cloud in the observable photosphere (r < 1); the lower value of implies the thicker cloud 
and hence a larger dust extinction, resulting in the smaller values of BCj. For the case B, 
the dust extinction is so large that the J flux is suppressed below the case of BB and hence 
the values of BCj tend to be smaller than those for the case of BB at Teflf ^ 1700 K. 




(4) 



Then 
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Next, we evaluate J — K based on the predicted fluxes of the UCMs, and the resulting 
J — K is reduced to the empirical system so that J — K = 0.0 for AOV star, again by 
the use of the model flux of Vega by Kurucz (1993). Wc used the filter response function 
of Persson et al. (1998) and their transform equation to the CIT system. The resulting 
Teff — (J — K) calibrations are shown in Fig. lb for the four values of Tcr, the cases B and C, 
and the case of BB. The red limit of J — K depends on Tcr and is larger (redder) for the lower 
value of Tcr, which implies the larger dust column density of the cloud in the photosphere; 
(J - K)^^ ^ 2.3, 2.0, 1.7 & 1.4 for T„ = 1700, 1800, 1850, & 1900 K, respectively (Fig. lb). 
On the other hand, the observed red hmit of J — X is about 1.9 from Fig. 4 of Dahn et al. 
(2002). Thus we suggest that ~ 1800 K in agreement with the previous analyses based on 
the different color systems (Tsuji 2001, 2002). This result is based on the grid of UCMs with 
log^r = 5.0, the changes of which by Alogf^ = ±0.5 result in A (J — K) ^ =1=0.2 (Fig. lb), 
and this may partly explain the scatter in the observed colors (Fig. 4 of Dahn et al. 2002). 

We now transform Mboi and T^s of the evolutionary models to Mj and J — K via Figs, la 
and lb with = 1800 K throughout. We first apply this procedure to the cooling tracks of 
the brown dwarfs with 10, 42, and 70Mjupiter by Burrows et al. (1997), and the results are 
shown in Fig. 2a with the observed Mj and J — K values (Dahn et al. 2002). It is to be noted 
that the red limit of L dwarfs as well as the rapid bluing of T dwarfs is well reproduced by the 
coohng tracks of the substellar mass objects. The position of the T2 dwarf SDSS J1254-0122 
is well accounted for by the cooling track of the moderate mass brown dwarf and this fact 
confirms that this object is in fact a transition object from L to T dwarf (Leggett et al. 2000). 
The high luminosity of the T5 dwarf 2MASS J0559-1404 at the J band is also accounted for 
by the cooling track of the low mass brown dwarfs. A few objects near Mj 15 may be 
accounted for by the models of higher masses while the observed photometric data of the 
very faint G1570D may not be well fixed yet. 

Next, we repeat the same analysis on the theoretical isochrones for 0.1, 1 and lOGyr 
by Chabrier et al. (2000) and the results are shown in Fig. 2b with the observed data. The 
observed position of 2M0559 is still higher than the isochrone of 0.1 Gyr, but a possibility 
of a binary was suggested (Burgasser 2002) and then the observed position may be shifted 
downward by 0.75 mag. The results shown in Figs. 2a and 2b are not fully self-consistent 
in that the model photospheres used in computing the evolutionary models are not the 
same as those used in converting the resulting (Mboi,Teff) diagrams to the CM diagrams. 
According to Burrows et al. (1997), however, by replacing gray models with nongray models, 
the evolutionary tracks have changed, but generally by no more than 10% in luminosity at 
any time, for any given mass. Also, the effect of the photospheric models on T^s was shown 
to be no larger than about 100 K (Fig. 2 of Chabrier et al. 2000). With these reservations 
in mind, we conclude that the general trend of the observed CM diagram can be reasonably 
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well fitted with the predicted ones based on the evolutionary models and the photospheric 
models (and Figs, la & lb can also be apphed to any other evolutionary track). 

It is to be noted that the "brightening" of the J fiux in the early T dwarfs is an 
artifact of observing objects with different masses and/or ages, and it should not imply that 
a single coohng track of a given mass shows the "brightening" . Nevertheless the rather high 
luminosity at the J band reflects the re-increase of BCj at the L/T transition noted before 
and this is a natural consequence of the migration of the dust cloud from the optically thin 
region (r < 1) in L dwarfs to the optically thick region (r > 1) in T dwarfs. This migration 
of the dust cloud recovers the gaseous opacities including H2 CIA at the K band region 
and makes the J flux large enough to compensate for the decreasing bolometric flux in the 
early T dwarfs. Also, the rapid transition from L to T dwarfs on the CM diagram is simply 
because the cloud moves from the optically thin to thick regions for only a small change of 
Teflf during the L/T transition. It is to be emphasized here that the migration of the dust 
cloud to the deeper photosphere in the cooler brown dwarfs is simply a consequence of the 
dust formation at the condensation temperature (Tcond ~ 2000 K) irrespective of T^s, and 
no other mechanism is needed at all to explain the migration of the dust cloud. 



3. DISCUSSION AND CONCLUDING REMARKS 

We have shown that the observed CM diagram extending from L to T dwarfs can be well 
fitted with the theoretical (Mboi, Tcs) diagram converted to (Mj, J — K) diagram via UCMs. 
This is the first successful application of the model photospheres to the CM diagram of the 
ultracool dwarfs which show rapid bluing and unexpected "brightening" at the J magnitude 
in T dwarfs after passing the red limit in L dwarfs. So far, other model photospheres 
including the cloudy models by Marley et al. (2002) as well as the dusty and clear models by 
Chabrier et al. (2000) could not be fitted to the observed CM diagram (Fig. 1 of Burgasser 
et al. 2002b). Although the cloudy models of Marley et al. (2002) have essentially the same 
feature as our UCMs in that the thin dust cloud is formed deep in the photosphere, yet their 
models failed to account for the observed (Mj, J — K) diagram, even if their free parameter 
referred to as the sedimentation efficiency, /rain, was changed. 

The difficulty of the Marley et al's cloudy models led Burgasser et al. (2002b) to suggest 
that the dust cloud should subject to disruption. Then the breaks between the clouds could 
explain the brightening of Mj as well as the bluing of J — of the T5 dwarf 2M0559, but 
only by introducing an additional free parameter, the cloud coverage fraction, besides the 
sedimentation efficiency, /rain- Also, the position of the L/T transition object SDSS 1254 
on the CM diagram could be explained only if the cloud coverage of 40% was assumed. 
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However, the argument based on the failure of the Marley et al's cloudy models to explain 
the observed CM diagram cannot be the evidence for the dust cloud disruption, since our 
simpler cloudy models, the UCMs, easily explain the observed (Mj, J—K) diagram including 
the key objects SDSS 1254 and 2M 0599 without any additional ad-hoc assumption, as shown 
in Sect. 2. As to an additional argument for the cloud disruption, namely the non-monotonic 
behavior of the FeH F^A-X^A (0 - 0) band at 0.9896 //m (Burgasser et al. 2002b), further 
observational confirmation will be needed since no FeH (1 - 0) band at 0.8692 /im can be 
seen in the spectrum of the T5 dwarf 2M0559 (Fig. 6.4 of Burgasser 2002). 

In conclusion, the observed CM diagram based on the recent parallax measurements 
can be fitted with the predicted CM diagrams based on the evolutionary models and UCMs 
of the photospheres. Since the theoretical evolutionary models and observed CM diagram 
can now be deemed as reasonably well established, this success in fitting the theoretical and 
observed CM diagrams can be regarded as an observational confirmation of the model pho- 
tospheres we have applied, namely the UCMs. This success of the UCMs may be significant 
especially because the UCMs are based on the very simple thermodynamical argument that 
the dust formed at T ^ Tcond can survive only to where T ^ Tcr, resulting in the formation 
of a thin dust cloud between T ^ T^^ and Tcond- Then, possibly complicated processes in- 
cluding condensation, growth, segregation, and precipitation of dust grains in the cloud are 
all absorbed in the Tcr, which is determined empirically. Such a semi-empirical approach 
plays an important role in stellar modeling, and a well known example of a useful empirical 
parameter is the so-called mixing length in the theory of stellar convection. Now our sim- 
plified treatment of the dust cloud has been shown to be successful for the interpretation 
of the CM diagram and will hopefully be useful for the interpretation and analysis of other 
observed data as well. 

We thank the anonymous referee for a careful reading and helpful comments. This work 
was supported by the grants-in-aid of JSPS Nos.11640227 (T.T.) and 14520232 (T.N.). 
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Fig. 1. — a) The values of BCj based on the unified cloudy models (UCMs) of log^r = 5.0 
with = 1700,1800,1850, and 1900 K are plotted against T^s by the solid lines. Also, 
those for the fully dusty models (case B), dust segregated models (case C), and blackbody 
radiation are shown by the dashed, dotted, and dash-dotted lines, respectively. The effects 
of changing \ogg by ±0.5 for the case of Teflf = 1600 K and = 1800 K are shown by the 
vectors, b) The same as for a) but for the values oi J — K. Note that the dust cloud results 
in essentially the same effect on BCj and on J — K (apart from the sign). 
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Fig. 2. — a) The cooling tracks of brown dwarfs with M = 10, 42, and 70 Mjupiter by Burrows 
et al. (1997) converted to CM diagrams via unified cloudy models (UCMs) of Tcr = 1800 K 
(solid lines) are compared with the observed data (filled circles) by Dahn et al. (2002). Also, 
shown are those with M = 42Mjupiter converted via the fully dusty models of case B (dashed 
line) and via the dust segregated models of case C (dotted lines), b) The same as for a) but 
for the isochrones of brown dwarfs for 0.1, 1.0 and 10 Gyr by Chabrier et al. (2000) converted 
to CM diagrams via UCMs of = 1800 K (solid lines). Also, those for 1 Gyr isochrone 
converted via the models of cases B (dashed line) and C (dotted lines) are shown. 



